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ABSTRACT 

Most models of anomalous X-ray pulsars (AXPs) account for the observed X-ray spectra and pulsations 
by means of radiation processes that occur on the surfaces of neutron stars. For any such model, general 
relativistic deflection of light severely suppresses the amplitude of the observed pulsations. We calculate 
the expected pulsation amplitudes of AXPs according to various models and compare the results with 
observations. We show that the high 70%) pulse amplitudes observed in some AXPs can be accounted 
for only if the surface emission is localized (spot radius < 40°) and strongly beamed cos" 9' with 
n ^ 2, where 9' is the angle to the normal). These constraints are incompatible with those cooling 
and magnetar models in which the observed X-rays originate as thermal emission from the neutron-star 
surface. Accretion models, on the other hand, are compatible with observations for a wide range of 
parameters. Finally, definitive conclusions cannot be reached on magnetospheric models, since their 
localization and beaming properties are not well understood. 

Subject headings: accretion, accretion disks — relativity — stars: neutron — X-rays: stars 

Submitted to The Astrophysical Journal. 



1. INTRODUCTION 

Among pulsating compact X-ray sources are a small sub- 
set with pulsation periods between 6 and 12 seconds, soft 
spectra, and no identifiable companions. The first sug- 
gestion that these objects might form a separate class of 
neutron stars - later called the anomalous X-ray pulsars 
(AXPs) - was made by Mereghetti & Stella (1995), who 
proposed that they might be powered by accretion from a 
very low mass companion (this was also earlier hinted at 
by Hellier 1994). The lack of optical counterparts, how- 
ever, and the absence of observable Doppler shifts in the 
frequency of the X-ray pulses led van Paradijs, Taam, & 
van den Heuvel (1995) to favor a different accretion model 
for AXPs. According to their suggestion, material from 
a fossil accretion disk, possibly the debris of a disrupted 
binary companion after a period of common-envelope evo- 
lution, is being accreted by a solitary neutron star. Re- 
cently, a similar model has been proposed by Chatterjee 
et al. (2000), in which the accreting material is supplied 
by the post-supernova fallback material from the neutron 
star progenitor itself. 

In a different class of models, AXPs are considered to be 
isolated neutron stars, spinning down by magnetic dipole 
radiation. Because of their unusually high period deriva- 
tives, a simple application of the relationship between the 
dipole magnetic field strength and the spin-down rate im- 
plies a field strength for these objects of ~ 10^^ — 10^^ G 
(see, e.g., Thompson & Duncan 1996). Two main types of 
models that rely upon the presumed high magnetic field 
of the stars have been proposed. Thompson & Duncan 
(1996; see also Duncan & Thompson 1992) suggested that 
the released energy may be drawn from the decay of the 
magnetic field itself and from differential movements in the 
stellar crust. This model also serves to explain the bursts 
observed from soft gamma-ray repeaters as being produced 



by larger-scale magnetospheric phenomena. Because of 
such models, AXPs are often called "magnetars" and are 
grouped into the same class of sources as the soft gamma- 
ray repeaters (see, e.g, Thompson & Duncan 1995, 1996; 
also Hurley 2000 for a review of SGR properties and mod- 
els). In the alternative high magnetic-field model, Heyl & 
Hernquist (1997a, 1997b) suggested that AXPs draw their 
energy from the residual thermal energy of the star itself. 

All these models face a number of difficulties. For exam- 
ple, if AXPs are powered by accretion from a stellar com- 
panion, the absence of detectable Doppler shifts in the ar- 
rival of X-ray pulses cannot be easily explained(Mereghetti 
et al. 1998). If accretion is from either a companion or a 
fossil disk, the optical fluxes one would expect directly 
from the disk or due to X-ray reprocessing are too high 
compared to the observed upper limits (see, e.g., Perna, 
Hernquist, & Narayan 2000; Hulleman et al. 2000). On 
the other hand, in the magnetar model the absence of 
bright faster AXPs and the observed variations in spin- 
down rates are hard to account for (see, e.g., Baykal & 
Swank 1996; Marsden et al. 2000; see, however, Heyl & 
Hernquist 1999; Melatos 1999). 

In all the above models of AXP, nearly all of the X-ray 
emission is produced at the surface of the neutron star. 
It is well known that strong gravitational flelds tend to 
smooth out the variability produced by a spinning com- 
pact star, even if the emission is highly localized in bright 
spots (Pechenick, Ftaclas, & Cohen 1983). Indeed, the X- 
ray pulse amplitudes of the three radio pulsars that show 
thermal emission from their surfaces are only ^ 30% (see, 
e.g., the discussion in Page 1995; Harding & Muslimov 
1998). This is in contrast to the non-thermal emission 
from radio pulsars (which is magnetospheric in origin) and 
from accretion-powered X-ray pulsars (which is from col- 
limated accretion columns) which often show pulse ampli- 
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tudcs as high as ^ 90% (see, e.g., Nagase 1989). In this 
respect, AXPs arc similar to accretion-powered pulsars, 
showing X-ray pulse amplitudes anywhere between ~ 10% 
and ~ 70%. 

Wc study in this paper a set of variability diagnos- 
tics that may be used in constraining emission models of 
AXPs. We examine three; parameters: the pulse fraction 
observed at infinity, which is a measure of the overall am- 
plitude of the variations, and the Fourier amplitudes at 
the first and second harmonics of the neutron-star spin 
frequency. We find tight constraints on the properties of 
magnetar models. Most importantly, we are unable to re- 
produce the observed variability properties of AXPs with 
thermal cooling models. 

2. FORMALISM 

In order to determine how relativistic effects suppress or 
enhance variability amplitudes, we need to consider curved 
photon paths from the surface of the star to an observer 
at infinity. Since the objects under consideration are ro- 
tating slowly, we use the Schwarzschild spacctime, which 
is appropriate for a non-spinning mass, and ignore effects 
such as relativistic frame dragging, which are important 
only for rapidly spinning objects. In this section, we out- 
line the basic ingredients of our method, drawing on the 
work of Pechenick et al. (1983). 

For each model we specify both the brightness distribu- 
tion over the surface of the star and the effective beam- 
ing of radiation, i.e., we specify the specific intensity inte- 
grated over photon energy, 1(0, (j), 9'), that emerges from 
each point on the stellar surface with polar coordinates 
{9, (j)) at an angle 9' with respect to the normal. 

In the current analysis we consider a number of differ- 
ent, physically motivated, mathematical expressions for 
the dependence of the emerging specific intensity on 9' 
(hereafter called the beaming function). This allows us to 
explore a large parameter space and draw conclusions that 
do not depend strongly on any particular emission model. 
We consider isotropic emission, i.e., no dependence on 9\ 
as well as the beaming described by the Hopf function 
(Chandrasekhar 1950, eq. [III.50]) 

/(m,^') = /o(m)(Ett^--^'+«) ' 

where the parameters La, ka, and Q are given in Chan- 
drasekhar (1950; Table III.VII). The Hopf function de- 
scribes the beaming of radiation emerging from a scatter- 
ing atmosphere heated from below and is, therefore, suit- 
able for a weakly-magnetic H- He atmosphere at energies 
^ 1 keV (Zavlin et al. 1998). 

For photon energies near the cyclotron energy {Ecyc — 
11.6[B/10^^ G] keV), the beaming function for a magnetic 
atmosphere does not always decrease monotonically away 
from the radial direction. Indeed the beaming function 
may have a local minimum at small angles from the nor- 
mal (Zavlin et al. 1995). However, for a dipole magnetic 
field this effect is not very significant and results in a rather 
small pulse fraction (< 30%) even when general relativistic 
effects are not taken into account (see Zavlin et al. 1995). 
Moreover, the cooling models of AXPs discussed here 
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Fig. 1. — Neutron star radii, in units of the Schwarzchild ra- 
dius, for different neutron star masses and equations of state (A: 
Pandharipande 1971; FPS: Lorenz, Ravenhall, & Pethick 1983; UU: 
Wiringa, Fiks, & Frabrocini 1988; L: Pandharipande & Smith 1975b; 
M: Pandharipande & Smith 1975a). 

require that the magnetic field of the neutron star is 
~ 10^^ — 10^"'' G and hence that the cyclotron energy is 
~ 1 — 10 MeV, i.e, much larger than the ~ 1 keV pho- 
ton energies that are of interest here. For such energies, it 
is reasonable to consider beaming fimctions that decrease 
monotonically away from the normal to the surface (Zavlin 
et al. 1995). 

In models of accretion columns, the interaction of ra- 
diation with a plasma in the strong magnetic field of the 
neutron star, as well as the possibility of radiation being 
obscured by the accretion column (see, e.g., Riffert et al. 
1993) leads to a much sharper beaming. For this reason, 
we also consider beaming functions of the form 

I{9,(t>,9') = Io{9,(t>)cos^9' . (2) 

According to Nagel (1981; see also Meszaros & Nagel 
1985), the radiation pattern emerging from an accretion 
column at low accretion rates is described quite well with 
n~2-3. 

Given a model for the X-ray emission from the stellar 
surface, the flux measured by an observer at distance d, 
whose polar coordinates with respect to the stellar rota- 
tion axis are (/?, $), is given by (cf. Pechenick et al. 1983) 

X / I{e,4>,9')xdxdy . (3) 

Jx=a Jy=0 

In the above equation 9, (j), and 6' depend implicitly on 
the angles x and y as described in Pechenick et al. (1983), 
Xmax = {R/M){1 - 2M/R)-'^/'^, M and R are the neutron 
star mass and radius, and we have set c = G = 1. The 
double integral has an integrable pole at Xmax- In the cal- 
culations presented here we have evaluated this integral to 
an accuracy of 10~^ using Romberg integration of the fifth 
order. 
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AXP PROPERTIES 
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References. — A: Oosterbroek et al. 1998 {SAX); B: Gottholf ct al. 
1999 {ASCA+SAX); C: Torri et al. 1998 {ASCA); D: Sugizaki et al. 
1997 {ASCA); E: Israel et al. 1999 {SAX); F: Iwasawa et al. 1992 
{GINGA). 

The degree of suppression of the pulsation ampHtude 
depends sensitively on the compactness of the neutron 
star. Figure 1 shows the ratio p = Rc^/2GM for different 
neutron-star masses and equations of state. Based on this 
figure, we limit our parameter study to p = 2, 3, and 4; 
larger values of p would correspond to unrealistically light 
neutron stars, even for the stiffest proposed equations of 
state (i.e., <^ 1.4M0 even for equation of state M). 
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3. PULSATION AMPLITUDES 

The pulsation amplitudes observed in anomalous X-ray 
pulsars allow us to place constraints on models of X-ray 
emission from their surface. In this paper, we do not at- 
tempt to fit particular observed pulse profiles but rather 
try to set general constraints on large classes of models. 
For this reason, we only consider the pulse fraction, de- 
fined as (cf. eq.B) 



PF 



rpmax. ZTimm 
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OO ~^ OO 



(4) 



and the Fourier amplitudes of the harmonics of the pulse 
frequency, defined by 



Foo (13.. $) = loiP) + hiP) cos($) + /2(/3) cos(2$) 



(5) 



Table I shows the observed pulse fractions and harmonic 
content for six AXPs, as inferred approximately from pub- 
lished lightcurves. We see that half of the known systems 
have pulse fractions of 0.5 or larger, with the largest be- 
ing ~ 0.7. As we discuss below, this fact provides strong 
constraints on some models of AXPs. 

3.1. Cooling of Magnetic Neutron Stars 

The brightness distribution on the surface of a strongly- 
magnetic cooling neutron star depends on the local mag- 
netic field strength B and its angle ip with respect to the 
local radial direction (see, e.g., Heyl & Hernquist 1998). 
The flux emerging from a spot on the stellar surface is 

i?™ cos^ ip, with m ~ 0.4 approximating well the numer- 



ical results (Heyl & Hernquist 1998). 
field we, therefore, use 

^ ,„ ^ cos^ 6 



(3cos2 6'+ l)i-™/2 



For a dipole stellar 



0.4, (6) 



where {0,(l>) are polar coordinates on the stellar surface 
with respect to the magnetic axis. 



Fig. 2. — Pulse fractions calculated for cooling neutron stars with 
a dipole magnetic field, and different magnetic inclinations (a) and 
orientations of the observer (/3) with respect to the rotation axis. 
The different panels correspond to an isotropic beaming function 
(n = 0), a Hopf function, and beaming functions that are increas- 
ingly more peaked towards the radial direction (n = 2, 3). 

Figure 2 shows the predicted pulse fractions for a cool- 
ing neutron-star model, for different choices of the angle a 
between the magnetic dipole axis and the rotation axis and 
the angle (3 between the light of sight to the observer and 
the rotation axis, as well as for different beaming functions. 
Some general features of this figure are worth noting. For 
example, the predicted pulse fraction remains unchanged 
when the magnetic inclination and the inclination of the 
observer are interchanged. This is true because the flux at 
infinity measured at any give pulse phase depends only 
on the angular distance 



= cos ^ (sin a sin /3 cos $ + cos a cos /3) 



(7) 



between the magnetic axis and the direction to the ob- 
server, which is symmetric in a and /?. Furthermore, the 
predicted pulse fraction is maximum when a — p. The 
overall maximum occurs when a — P — 90° and the pulse 
fraction is zero when either a — 0° or /3 — 0°. 

The most important result from Figure 2 is that even 
for the most favorable geometry (i.e., a = P = 90°) and 
for strong beaming (i.e., n = 3), the pulse fraction does 
not exceed 37%. The reason is that the surface brightness 
distribution (o) is too smooth. Assuming a less relativistic 
neutron star (i?/2M = 4), or even a significantly stronger 
dependence of the emerging flux on the local magnetic field 
(m = 1 in eq. ^) results in only a modest increase of the 
maximum pulse fraction (51% and 50%, respectively, for 
n = 3). We thus conclude that neutron star cooling mod- 
els with the surface emission described by equation (|^) 
cannot reproduce the large pulse fractions observed from 
AXPs. 
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Fig. 3. — Ratio of the bolometric flux emerging from an atmo- 
sphere with a given mass fraction {AM /M){Q/ An) (in terms of the 
neutron-star mass M) of H-He matter accreted over a fraction Q/4it 
of the steUar surface, to the bolometric flux from a purely iron atmo- 
sphere with the same core temperature (after Chabrier et al. 1997). 



3.2. Localized Thermal Emission 

Thermal emission from the neutron-star surface can, 
in principle, be more localized than indicated in equa- 
tion , if it is confined mostly around the stellar magnetic 
poles. For example, metallicity gradients on the surface of 
a cooling neutron star, possibly produced by magnetically- 
channelled fallback material during the supernova explo- 
sion, can lead to larger effective temperatures near the 
magnetic poles than in the magnetic equator (see, e.g., 
Pavlov et al. 2000) . Alternatively, non- uniform heating of 
the neutron-star atmosphere, e.g., caused by magnetic field 
decay or crustal heating (see, e.g., Thompson & Duncan 
1996), can lead to more localized thermal emission from 
its surface. 

We estimate the magnitude of the first effect using the 
analytic expressions for neutron-star atmosphere models 
given by Chabrier, Potekhin, & Yakovlev (1997). We as- 
sume that a fraction (AM/Af)(il/47r) of light element ma- 
terial has accumulated only over a fraction (i^/An) of the 
stellar surface. We neglect the fact that the surface lay- 
ers of the neutron star are in the liquid phase and hence 
lateral diffusion may smooth out the metallicity gradient 
(L. Hernquist, private communication). For a given core 
temperature we then calculate the ratio / of the bolo- 
metric flux emerging from the light-element region of the 
surface to the bolometric flux from the region consist- 
ing purely of iron. This flux contrast is shown in Fig- 
ure 3, for different values of the mass fraction of the ac- 
creted atmosphere. The flux contrast does not further 
increase, for (AM/M)(r2/47r) > 10^^, because the base of 
the accreted layer reaches densities that are high enough 
(~ 10^° g cm~^) for it to be part of the isothermal core. Al- 
though the above models are strictly valid only for weakly 
magnetic neutron stars, the flux ratios are significantly 
more sensitive to the composition of the neutron-star en- 
velope than the magnetic field strength (cf. Heyl & Hern- 
quist 1997b). 
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Fig. 4. — Calculated pulse fractions for a neutron star with 
R/2M = 3, as a function of the half-opening angle of each polar 
cap, for different beaming functions, and a flux contrast of 9 be- 
tween the polar caps and the rest of the stellar surface. 

According to Figure 3, the expected flux contrast be- 
tween the polar caps and the rest of the stellar surface can 
be at most / ~ 9 and the contrast depends very weakly 
on the amount of accreted material. For the flux contrast 
to attain its maximum value, the polar caps must accrete 
AM ^ 10~^Mq, while the rest of the surface must accu- 
mulate more than three orders of magnitude less material. 
Even if these conditions were met and even if the photons 
emerged from the stellar surface strictly radially, the pulse 
fraction would be at most (/ — l)/(/ + 1) 0.8. 

Any realistic beaming function reduces this upper bound 
significantly below the highest observed pulse fraction for 
AXPs (~ 0.7; see Table I), even for mildly relativistic neu- 
tron stars. This is demonstrated in Figure 4, which shows 
the maximum pulse fraction of radiation emerging from a 
neutron star with R/M = 3, as a function of the angular 
radius of each polar cap, for various realistic beaming func- 
tions, and an emerging flux from the polar caps that is nine 
times larger than the rest of the stellar surface. For small 
polar-cap sizes, the radiation flux emerging from the caps 
is only a small perturbation to the total brightness of the 
star and hence the pulse fraction is small. With increasing 
cap size the relative contribution of the caps to the bright- 
ness of the star increases, leading to an increase of the 
pulse fraction, until the caps cover a large enough fraction 
of the stellar surface and the pulse fraction drops again. 
For the flux contrast and neutron star radius used here, 
the maximum pulse fraction occurs when the polar-cap 
size is ~ 40° and this maximum value is 55%. Clearly, 
no realistic model of this kind can fit the pulse fraction of 
~ 0.7 seen in IE 1048.1-5937. 

The pulse fraction can be further enhanced if the fluxes 
emerging from two antipodal polar caps are unequal, e.g., 
because of uneven fallback or crustal heating. As the 
contrast between the two polar caps increases, the pulse 
fraction becomes larger, but only at the expense of the 
Fourier amplitudes at the harmonics of even order. This 
is shown in Figure 5, where the predicted ratio I2/I1 is 
plotted against the pulse fraction for R/M = 3, two 
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Fig. 5. — Ratio l2/I\ of the Fourier amplitudes at the second 
harmonic and fundamental of the pulsar spin frequency (cf. eq. H) 
as a function of the total pulse fraction (cf. oq. M). The calculations 
correspond to R/M = 3, a configuration of two antipodal polar caps 
with a half-opening angle of 40 degrees but different emerging fluxes, 
and beaming described by the Hopf function. The solid lines outline 
the area in the parameter space allowed for specific ratios of fluxes 
emerging from the two polar caps. The dashed line represents the 
maximum ratio I2/I1 that corresponds to a given pulse fraction; for 
such a configuration, no system is allowed to lie to the right of the 
dashed line. The diamonds correspond to the observed properties 
of AXPs (cf. Table 1). 

antipodal caps with a half-opening angle of 40 degrees, and 
beaming described by the Hopf function. For these calcu- 
lations, the two caps are assumed to have different emerg- 
ing radiation fluxes, with the flux of the brightest set to 9 
times the flux emerging from the rest of the neutron-star 
surface. For any given pulse fraction, there exists a max- 
imum value that the ratio I2/I1 can attain, shown by the 
dashed line in Figure 5. As a result, detection of a source 
with a large pulse fraction and significant amplitude at the 
even order harmonics can exclude such a configuration. 

The observed properties of some AXPs are not consis- 
tent even with a model with such unequal polar caps (Fig- 
ure 5). For example, the source IE 2259-1-586 is character- 
ized by a large pulse fraction and high harmonic content 
that cannot be achieved by any of the configurations con- 
sidered here. It does not, therefore, appear plausible for 
thermal emission from the stellar surface to account for 
the variability properties of AXPs. Achieving the kind of 
variability observed is possible only when the emission is 
both localized and strongly beamed. 
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Fig. 6. — Calculated pulse fractions for a neutron star with 
R/2M = 3, a hot-spot half-opening angle of 5 degrees and a cos" d' 
beaming function, with n=0— 3. The maximum value of the pulse 
fraction is written above each panel. 
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3.3. Accretion Onto Magnetic Neutron Stars 

In accretion models of AXPs (see, e.g., Mereghetti & 
Stella 1995; van Paradijs et al. 1996; Chatterjee et al. 
1999), a large fraction of the X-ray emission may be pro- 
duced mainly at localized "hot spots" where the accre- 
tion columns meet the stellar surface. In order to take 
such configurations into account, we describe the surface 
brightness distribution with two circular antipodal spots, 
with a brightness that is constant over their surface area. 
We denote by a the angular distance of the center of each 
spot from its closest rotation pole. 



Fig. 7. — Calculated pulse fractions for a neutron star with 
R/2M = 3, a cos'^ 6' beaming function, and different hot-spot half- 
opening angles. The maximum value of the pulse fraction is written 
above each panel. 

If AXPs are powered by accretion from a geometrically 
thin disk, the half-opening angle of each polar cap can be 
very small. We estimate the size of each polar cap as the 
angular distance on the stellar surface from the magnetic 
axis of the footpoint of the last magnetic field line that 
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Fig. 9. — Contours of maximum pulse fraction for different beam- 
ing functions (cos" 9'; cf . eq. H] ) and half-opening angles of each hot 
spot. The ~ 70% pulse fraction observed from the source IE 1048.1- 
5937 restricts any viable model of AXPs to the left of the uppermost 
contour. 




Fig. 8. — Calculated pulse fractions for different neutron-star 
radii (p = R/2M), a hot-spot half-opening angle of 5 degrees, and a 
cos^ 8' beaming function. The maximum value of the pulse fraction 
is written above each panel. 

penetrates the accretion disk. Given that all AXPs are 
observed to be spinning down, their accretion disks must 
be truncated near the outer corotation radius (see, e.g., 
Ghosh & Lamb 1979) 

-Rco = I 
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:550i? 



/ R 



V 106 cm 



M 



1AM, 
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6^ 
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where P is the spin period of the pulsar. For a dipole mag- 
netic field, the quantity sin^0/r, with 9 measured from 
the magnetic pole, remains constant along a field line, 
and therefore the half-opening angle of the polar cap is 
^ {R/Rco)'/^ ~ 2.5°. 

Figure 6 shows the pulse fraction calculated for a model 
with two identical hot spots, a half-opening angle of 5 de- 
grees, and various beaming functions (see eq. Q). Because 
of the general relativistic deflection of light, we see that 
even such a small spot cannot produce a ^ 70% pulse frac- 
tion, unless there is significant beaming (n ^ 2). On the 
other hand, if the emerging radiation is strongly beamed 
towards the radial direction, a large pulse fraction can be 
achieved for a wide range of polar cap sizes (<; 40°; Fig. 7) 
and for all realistic neutron-star radii (Fig. 8). 

4. CONCLUSIONS 

In this paper wc have used the high pulse amplitudes ob- 
served from a number of AXPs to constrain the properties 
of their emission mechanism. We find that the observa- 
tions can be accounted for only if the surface emission is 
localized (half-opening angle < 40°) and strongly beamed 



(n ^ 2 in eq. Q), as summarized quantitatively in Fig- 
ure 9. These constraints are a consequence of the compact- 
ness of the neutron stars and the resulting strong general 
relativistic deflection of photon trajectories. Our conclu- 
sions are valid for all realistic neutron star masses and 
radii. 

The properties of individual sources offer a number of 
additional clues. For example, the double-peaked pulse 
profile of IE 2259-1-586 requires that the emission is local- 
ized around two antipodal spots on the neutron-star sur- 
face, probably associated with the magnetic poles. Fur- 
thermore, the change in the relative strength of the two 
peaks observed with GINGA (Iwasawa et al. 1992) implies 
that the pulse shape cannot be solely due to geometric ef- 
fects but should also reflect a flux contrast between the two 
antipodal spots. Furthermore, this flux contrast should be 
variable and, therefore, cannot be caused by a non-dipolar 
magnetic field configuration or a non-uniform fallback of 
low metallicity material. 

Such arguments, together with the constraints presented 
in Figures 2, 4, 5, and 9 appear to rule out thermal cooling 
models for AXPs. They are also inconsistent with those 
magnetar models in which most of the X-ray flux origi- 
nates from heating in the deep surface layers of the neu- 
tron star. On the other hand, the localized emission and 
beaming predicted by accretion models seem to be con- 
sistent with the observations. A magnetospheric model 
in which the neutron-star surface is heated over localized 
spots by particle bombardment may also be viable, though 
the beaming properties of such a model are unknown. 

We thank Deepto Chakrabarty, Lars Hernquist, Vicky 
Kaspi, Jessica Lackey, and Feryal Ozel for many useful dis- 
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doctoral fellowship of the Smithsonian Institution. 
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